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ABSTRACT 

The P/ancfc-ATCA Co-eval Observations (PACO) project collected data between 4.5 
and 40 GHz for 482 sources selected within the Australia Telescope 20 GHz (AT20G) 
catalogue and observed with the Australia Telescope Compact Array (ATCA). Ob- 
servations were done almost simultaneously with the Planck satellite, in the period 
between July 2009 and August 2010. In this paper we present and discuss the data for 
the complete sample of 159 sources with S'at20G > 200 mJy in the Southern ecliptic 
pole region. 

The Planck Early Release Compact Source Catalogue (ERCSC) contains 57 of our 
sources. A comparison between the PACO catalogue and the ERCSC confirms that 
the reliability of the latter is better than 95 per cent. The missing ERCSC sources are 
typically associated with the Large Magellanic Cloud, the Milky Way or are otherwise 
extended. 

The spectral analysis of the PACO faint catalogue shows a spectral steepening of 
the sources at high frequencies, confirming the results obtained from the PACO bright 
sample. 

A comparison with AT20G measurements, carried out, on average, a few years 
earlier, has demonstrated that, on these timescales, our sources show a rather high 
variability with an rms amplitude of ~ 40 per cent at 20 GHz. The source spectral 
properties are found not to vary substantially with flux density, except for an increase 
of the fraction of steep spectrum sources at fainter flux densities. 

Our data also allow us to extend by a factor ~ 5 downwards in flux density 
the source counts at ~ 33 GHz and ~ 40 GHz obtained from PlancKs Early Release 
Compact Source catalogue. This allows us to substantially improve our control on the 
contribution of unresolved extragalactic sources to the power spectrum of small scale 
fluctuations in CMB maps. 

Key words: galaxies: active - radio continuum: galaxies - radio continuum: general 
- cosmic microwave background. 



1 INTRODUCTION 



The knowledge of the extragalactic radio sources at 
frequencies higher than 10 GHz is still poor, al- 
though the situati on has been impr oving rapidly in 
recent years (see |de^ Zotti et all |2010|. for a review), 
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IOC (iDavies et all I2OI0I') WMAP (IWright et all l2009l : 
iLopez-Can icgo ct al."2007'. '2009'; 'Massardi et al.ll2009l ') , and 
Planck ([Planck Collaboration 2011a) surveys. The study of 
extragalactic sources at millimeter wavelengths is impor- 
tant both per se and in connection with CMB experiments, 
since these sources are the main astrophysical contaminant 
of CMB maps on angular scales small er than 20'-30' at fre- 
quencies of up to ~ 100 GHz (|Toffolatti et al.l Il998l . Il999l : 
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Planck offers a unique opportunity to carry out an unbi- 
ased investigation of the spectral properties of radio sources 
in a poorly explored frequency range, partially unaccessible 
from the ground. A definition of sources Spectral Energy 
Distributions (SEDs) over a frequency range as large as pos- 
sible is crucial in determining their physical properties and 
in identifying the different components that may contribute 
to their emission. Given that observations in the full Planck 
frequency range will not be repeatable at least in the fore- 
seeable future, it is essential not to lose the occasion, while 
Planck is flying, of greatly increasing its science yields by 
simultaneous (i.e. not affected by variability) ground based 
observations at lower frequencies as well as at frequencies 
overlapping with Planck channels. 

This has motivated the Planck-ATCA. Co-eval Obser- 
vations (PAGO) project, that consisted in observations with 
the Australia Telescope Compact Array (ATGA) at several 
frequencies between 4.5 and 40 GHz, almost simultaneously 
with the Planck satellite. In the period between July 2009 
and August 2010 we observed 482 extragalactic sources, ex- 
tracted from the Australia Telescope 20 GHz (AT20G) sur- 
vey catalog (Murphv ct al. 2010). Of these, 344 sources form 
3 partially overlapping complete sub-samples, selected for 
different purposes. 

• The PAGO faint sample, presented in this paper, is 
made of 159 sources with 52oghz ^ 200 mjy in the South 
Ecliptic Pole region (ecliptic latitude < —75°) and with 
3h < RA < 9h, 5 < -30°. Near the Ecliptic Poles Planck's 
scan circles intersect. Therefore the area is covered many 
times, and Planck's sensitivity is maximal. 

• The PAGO spectrally-selected sample comprises the 69 
sources with S20GHZ > 200 mJy and inverted or upturning 
spectra in the frequency range 5-20 GHz, selected over the 
whole Southern sky; 

• The PAGO bright sample comprising the 189 sources 
with 5'20GHz > 500 mJy at 5 < -30° (Massardi ct al. 2011). 

A full desc ription of the PAGO p roject and of its main goals 
is given in iMassardT et al.l (|201lh . The aims specific to the 
PAGO faint sample are: 

• Extend to fainter flux densities the characterization of 
radio source spectra from 4.5 GHz to the Planck frequency 
range. Although Planck's sensitivity is maximal in the se- 
lected area, we expect that only a minority of sources in 
this sample will be detected by Planck. The estimated com- 
pleteness limits o f the Planck Early R elease Point Source 
Gatalog (ERGSG. IPlanck Gollaboration 2011c) are of 1 and 



1.5 Jy at 33 and 40 GHz, respectively (iPlanck G ollaboration! 
bOllil ). and the minimum flux density of an ERGSG extra- 
galactic source is of 480 mJy at 30 G Hz and of 585 mJy at 
44 GHz l|Planck Gollaboration|[201ld ). Only 57 PAGO faint 
sources have ERGSG flux density measurements at least at 
one frequency. The final catalog is ex pected to reach deeper 
flux density limits (|Leach et al. |2008| ) but, nevertheless, we 
will probably get the full Planck spectral coverage only for 
a subset of sources in this sample. However further spectral 
information can be obtained via a stacking analysis. 

• Extend the determination of source counts at ~ 33 
and ~ 40 GHz obtained fro m the analysis of the ERGSG 
llPlanck Gollaborationll2011il ') downwards in fiux density by 
a factor of ~ 5. Going down in flux density is important 



to control the contamination of GMB maps by faint radio 
sources. 

The plan of the paper is as follows. In §[2] we summarize the 
main steps of data acquisition, reduction and calibration. In 
§|3]we describe the PAGO faint sample. In §[3]we analyze the 
spectral behaviour between 4.5 and 40 GHz of our sources, 
their variability, and estimate the differential source counts 
at 33 and 40 GHz. Our main results are summarized in §[3 



2 OBSERVATIONS AND DATA REDUCTION 

A detailed description of o bservations, data redu ction and 
calibration can be found in iMassardi et al] (|201ll ). Here we 
present only a brief summary of the main points. The obser- 
vations exploited the capabilities of the new Gompact Array 
Broadband Backend (GABB, Ferris & Wilson 2002, Wil- 
son et al. 2011) system, that allows to observe 2x2 GHz 
simultaneous bands in continuum. Each band is split into 
2048 X IMHz frequency channels. Applied to the 6 x 22 m 
antennas of the Australia Telescope Gompact Array, the 
GABB gives a noise level per band reaching 0.5 mJy in 1 
min on source down to the 7 mm wavelength band. We chose 
to use the receivers at 7mm (intermediate frequencies, IPs, 
at 33 and 39 GHz), to overlap the lower- frequency Planck 
channels, at 12mm (IPs at 18 and 24 GHz) to overlap our 
selection frequency, and at 3-6 cm (IFs at 5.5 and 9 GHz) 
to extend the SEDs to lower frequencies. 

Since most of the sources are point-like we can observe 
with any array configuration. The dates of the observing 
runs, the telescope configuratio ns and their properti es are 
summarized in Tables 1 and 2 of lMassardi et ahl (|201ll ) . The 
project was allocated 62 observing runs, for a total of ^ 450 
hours. For each run and for each frequency, we observed a 
bandpass calibrator and a primary calibrator. Thanks to the 
compactness of most of the PAGO sources neither imaging 
nor phase calibration is necessary, but since sources are self- 
calibrated we get good images when observing with hybrid 
arrays. We got a suitable flux density estimation from the 
visibilities using the triple-correlation techniques. Each tar- 
get source has been observed in a single 1.5 min- long shot 
for each frequency, that corresponds to a theoretical noise 
of less than 1 mJy. 

Observations were generally made within 10 days from 
the Planck scan on any given source. An effort was made to 
observe sources in all the three bands (7mm, 12mm and 
3-6cm) during the same run. Exceptions are due to bad 
weather conditions. The sources have been re-observed in 
several runs overlapping with the Planck observations (the 
satellite surveys the whole sky every 6 months), with a 
higher observing rate for the sources that appeared more 
variable. 

For data reduction we u sed the MIRIAD software 
(jSault. Teuben. fc Wright! ri995h . recently updated to deal 
with GABB data. We inscribed MIRIAD tasks within a 
scripting pipeline for consistency in the analysis, and have 
independently reduced the data from each observing run and 
each frequency. After flagging the time intervals affected by 
bad weather or instrumental failures we have split the ~2 
GHz bands in 4x512 MHz sub-bands to exploit the spectral 
information. We have calibrated each sub-band for band- 
pass, ffux density and leakages. 



The PACO faint sample 3 



The classification of sources as point-like or extended 
was mostly taken from the AT20G catalogue and integrated 
at 7 mm by the analysis of the source phase closure and ra- 
tio between the data collected with long baselines and those 
collected with the whole array. The flux densities, and the 
associated errors, for point-like sources were computed usin g 
the triple product (|Thompson. Moran. fc SwensonI l200lh . 
For extended sources we used the scalar flux density on the 
shortest baseline which, however, is only a lower limit to 
the real flux density. We have in total 9 extended sources in 
the PACO faint sample, only 4 of which were classified as 
extended in the AT20G catalog. 

At the end of the process, a quality check on the data 
has been performed comparing the data points with a poly- 
nomial fit to the source SEDs. This approach assumes that 
the SEDs of our sources are smooth across our frequency 
range. Misbehaving points are attributed to frequency de- 
pendent problems, not addressed by the automatic flagging 
that operates within each bands. 



3 THE SAMPLE 

The region covered by the PACO faint sample is shown 
in Fig. [1] The crosses are the PACO faint sample sources, 
the circles are the 30 GHz ERCSC sources. There are 
16 ERCSC sources within the region bounded by the 
dashed lines, not included in the present sample. Twelve 
lie either within the excluded Large Magellanic Cloud 
(PLCKERC030 G276.20-33.22, PLCKERC030 G277.16- 
36.06, PLCKERC030 G277.62-32.14, PLCKERC030 
G278.40-33.32, PLCKERC030 G279.49-31.67) or |fo| < 5° 
(PLCKERC030 G252.05-h03.28, PLCKERC030 G260.67- 
03.20, PLCKERC030 G263.73-03.72, PLCKERC030 
G265.14-h01.48, PLCKERC030 G267.93-01.02, PLCK- 
ERC030 G253. 63-00.23, PLCKERC030 G248.44-04.00) 
areas. Three (PLCKERC030 G271.21-09.56, PLCKERC030 
G249.07-05.17, PLCKERC030 G251. 04-05. 33) are at low 
Galactic latitude and were not detected by the AT20G 
survey, suggesting that they are extended Galactic sources 
or peaks of the Galactic diffuse emission. The last one 
(PLCKERC030 G240.02-56.83) is Fornax A, missed by 
the AT20G because its core emission at 20 GHz is below 
the detection limit and its extende d emission has a siz e 
exceeding the ATCA field-of-view l|Massardi et all l200l ). 
Thus the fraction of potentially spurious ERCSC sources in 
the PACO faint area is < 5%. 

The results of the cross-match between PACO and 
ERCSC sources in all Planck channels are summarized in 
Tab. [U where we also list the adopted search radii and the 
rms values of the distances between Planck and PACO po- 
sitions in each channels. The rms values are much lower 
than the search radii, and consistent with the rms scat- 
ter in positions being lower than FWHM/5 as found in 
iPlanck Collaboration! l|2011 j ). Fig. [2] shows the distribution 
of the angular distances between the positions of the sources 
obtained by cross-matching the ERCSC catalogue with the 
PACO faint sample one: we chose the 100 GHz (grey area) 
and 143 GHz (hatched area) Planck channels, where we 
found the greatest number of ERCSC-PACO common ob- 
jects. Of the 57 sources present in the ERCSC, 25 have 
measurements at more than 5 different frequencies. While 



o' oo 



oh :0 o "SoO Oq 

o ^ o ; 

-10 - \ ^ 



oi o 

O : O 



O 

e' 



o 



o 



.cT.i.x .0.:. ..igp ^OH 

X V Q 



..■a. 



o-OQ 



300 290 280 270 260 250 240 230 

Glon [deg] 

Figure 1. The area covered by the PACO faint sample is that 
bounded by the dashed Unes, excluding the grey areas, i.e. the 
Galactic plane region (|6| 5°) and the Large Magellanic Cloud 
(LMC) region (a circle of 4° radius). The crosses show the po- 
sitions of PACO sources, the circles those of 30 GHz ERCSC 
sources. 



Table 1. Number of PACO faint sources present in the ERCSC, 
for each frequency channel. The last column gives the rms dis- 
tances between ERCSC and PACO positions for the channels with 
a number of sources greater than 20. 



freq 
(GHz) 


search radius 
(arcmin) 


number of 
sources 


rms distance 
(arcmin) 


30 


15 


39 


2.6 


44 


12 


27 


2.1 


70 


7 


29 


1.3 


100 


5 


42 


1 


143 


4 


44 


0.8 


217 


3 


31 


0.3 


353 


3 


15 




545 


3 


5 




857 


3 


3 





a joint analysis of PACO and ERCSC results is deferred to 
a subsequent paper, we show in Fig. [3] some examples of 
spectra determined combining the two data sets. 

All the 159 sources in the PACO faint sample have 
been observed at least once at all the frequencies men- 
tioned in §[21 Bad weather conditions affect all the runs for 5 
sources (J034205-370320, J040848- 750720, J052006-454745, 
J062620-534151, J064118-355433) and it was not possible to 
repeat the observations. All of them are extended. More- 
over, the results of mapping and mosaicing in total intensity 
and polarisation the core and the western lobe of Pictor A, 
namely J051926-454554 and J051949-454643, will be pre- 
sented in a subsequent paper; these two sources will not 
be considered further here. Therefore the PACO faint cata- 
logue, which is available as online supplementary material, 
contains 152 sources. The content of its columns is as fol- 
lows. 

1. AT20G name; 

2. -3. J2000 equatorial coordinates; 
4. Epoch of observation; 
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Figure 3. Examples of radio spectra obtained combining PACO (asterislss) and ERCSC (diamonds) measurements. 
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Figure 2. Distribution of the angular distances (in arcmin) be- 
tween the PACO and ERCSC positions of the sources resulting in 
the cross-match with the 100 (grey area) and 143 GHz (hatched 
area) Planck channels. 



5. Flag 's' for simultaneity with Planck; 

6. Flag 'e' for extended sources; 

7. -31. Flux densities in mjy for the 24 PACO sub-bands; 



32.-55. Flux density errors; 

56.-60. Double power-law fit parameters (see §|31). 



4 DATA ANALYSIS 

The analysis was restricted to the 143 point-like sources. The 
PACO data (4.5-40 GHz) have been fitted with the double 
power-law 



S{v) 



(1) 



with V and vq in GHz, S and So in Jy. S(i,vo,a and b are free 
parameters. For spectral indices we adopt the convention: 



log(5'l/g2) 



(2) 



The fit was performed for sources having more than 4 data 
points for each receiver band. This was the case for all but 2 
sources. The fit was considere d successful if the reduced x < 
2.9 (see lMassardi et aLllioTH ). The fit was not successful for 
nine sources which show an irregular spectrum that could 
not be described by the fitting formula in eq. ([l|. In most 
cases the irregularities of the spectrum seem to be due to 
residual calibration problems. When observations at more 
than one epoch are available, we chose the one with the 
greatest number of data points. 
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Figure 4. Distribution of PACO faint sources with 200 mjy ^ 
'Sat20G < 500 mJy in the Og''— alj] plane. The dotted lines show 
the boundaries for the various spectral types. 



Table 2. Percentages of sources of each spectral type in the 
PACO faint sample (200 Sat20G < 500 mJy) and in the bright 
sample (5at20G ^ 500 mJy). 



Type 200mJy ^ S'20GHz < 500mjy 
per cent 



5'20GHz 500mJy 
per cent 



flat 


5.1 


10.3 


steep 


13.3 


3.6 


inverted 





0.6 


peaked 


11.2 


14.5 


down turning 


65.3 


66 


self absorbed 


5.1 


4.8 


upturning 









4.1 Spectral behaviour 

The double-power-law fit that we have performed allows 
us to estimate the spectral indices for any frequency range 
between 4.5 and 40 GHz. Excluding sources brighter than 
500 mJy, which are included in the PACO bright sample, 
whose spectral properties are discussed by iMassardi et al.l 
we have 98 sources in the range 200 mJy ^ S'at20G < 
500 mJy. The spectral classification, based on the spectral 
indices between 5 and 10 GHz, compared with those be- 
tween 30 and 40 GHz {af , agg), is illustrated in Fig. g] 
The fractions of sources of each spectral type are listed in 
Table [2l where the fraction found for the bright PACO sam- 
ple are also given, for comparison. The main difference with 
the bright sample is the larger fraction of steep-spectrum 
sources, mostly at the expense of peaked- and flat-spectrum 
sources. 

For the 11 peaked sources with 200 ^ 5'2oghz < 
500 mJy we computed the peak frequency as Vp — 
uo{~h/aY^^^^-''^\ we found a median value of 16.6 GHz, 
a mean of 17.2 GHz and a dispersion of 3.7 GHz. For the 
64 down-turning sources in the same flux density range we 
found a median break frequency value of 31.8 GHz, a mean 



Figure 5. Distributions of Og" (shaded area) and of (hatched 
area). 
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Figure 6. Distribution of the difference ctg" — o|q. 



of 26.2 GHz and a dispersion of 14.1. We note that 28 sources 
have a value of vq at the border of our frequency range and 
therefore not reliably estimated. This does not affect the 
value of the median. 

The distributions of spectral indices Qg" and are 
shown in Fig. [5] The median values are -0.04 for af and 
—0.69 for Q30, with interquartile ranges [—0.22; 0.15] and 
[—0.93; —0.43] respectively, showing that the spectra are 
mainly flat at lower frequencies and steepen at higher fre- 
quencies. The distribution of the differences between the 
two spectral indices shown in Fig. [6] has mean value of 
0.67 ± 0.04 and standard deviation of 0.35; the median 
value is 0.63 ± 0.04. A steepening of the spectr al indices be- 
tween 20 and 90 GHz has also been reported bv lSaiina et al.l 
based on VLA observations of a sample selected from 
the AT20G survey. Moreover a steepening of spectral in- 
dices at higher freque ncies is also shown in the a nalysis of 
the NEWPS sar nple llGonzalez-Nuevo et al. l liooil ') and the 
ERCSC sample (|Planck Collaborationll2oTTl ). 

We did the same analysis for the spectral indices be- 
tween 5 and 20 GHz and between 20 and 40 GHz. The mean 
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Figure 7. Comparison between the AT20G flux densities at 20 
GHz and the PACO ones scaled to 20 GHz. The solid line corre- 
sponds to S'at20G = S'pACO a-nd the dot-dashed lines correspond 
to SpACO = SaT20G,20GHz(1 ± <^v)- 



value of Qg is —0.16 with standard deviation of 0.31; the 
mean value of q|q is —0.60 with a standard deviation of 0.31. 
These values are similar, although not equal to those found 
for the bright sample (—0.07 and —0.55, respectively, with 
standard deviations of 0.32 and 0.34). This indicates that 
the overall spectral behaviour up to 40 GHz is not strongly 
dependent on flux density in the range covered by these 
samples, especially between 20 and 40 GHz. 

4.2 Variability 

The AT20G flux densities have been collected in several 
epochs between 2004 and 2008. Therefore, the time lag be- 
tween the AT20G and the PAGO observations ranges from 
several months to several years. Since variable sources are 
more likely to meet the detection limit of a blind survey if 
they are in a bright phase, it is not surprising that, on aver- 
age, AT20G flux densities are slightly higher than the PACO 
ones. After correcting for the small frequency difference be- 
tween the PAGO and AT20G observations, using our mea- 
sured spectral indices between 18.768 23.232 GHz, we flnd 
(logio(S'AT20G/S'pAco)> = 0.043, where 5at2og and Spaco 
are the AT20G flux density at 20 GHz and the PAGO flux 
density scaled to 20 GHz. A comparison of the two flux den- 
sities is shown in Fig.[71 The least-squares linear relationship 
is SpAco[Jy] = (0.89 ± 0.04) ■ SAT20G[Jy] + (0.05 ± 0.05), 
showing no statistically signiflcant deviations from a con- 
stant ratio between the two flux densities. It is important to 
note, in this respect, that the observations have been per- 
formed with the same front-end system, but the back-end 
has changed. Last, but not least, the new GABB system 
requires a more careful account of the spectral behaviour 
of the sources during the analysis, because of the broader 
bands. Our flndings add to a series of indications of the 
high performance of the new system and show that there 
are no obvious instrumental effects that may complicate the 
comparison with measurements obtained before its upgrade. 

Since both the PAGO and the AT20G measurement er- 
rors are negligibly small, the variability at ~ 20 GHz can 




Figure 8. Comparison between the AT20G flux densities at 4.8 
GHz and the PACO ones at 4.744 GHz. The solid line corre- 
sponds to 5at20G = SpACO E'lid the dot-dashed lines correspond 
to S'paCO,4.7GHz='S'aT20G,4.8GHz(1 ± 0"d)- 



be estimated as the standard deviation of the distribution 
of S'paco/5'at20G or of logio(5'pACo/5'AT20G)- The distri- 
bution of the latter quantity is better fitted by a Gaussian 
than that of the former one, see Fig [9] The dotted line is 
the Gaussian best-fit, centered at -0.03 with FWHM = 0.34, 
and the dashed lines are the mean and the mean ± cr of the 
distribution. The mean and the standard deviation, o, of 
log^Q(5'pACo/'S'AT20G) are —0.043 and 0.18, respectively, cor- 
responding to a dispersion (t„ — cr (5'paco/'S'at2og) — 0.40, 
on a timescale of a few years. For comparison, for the bright 
sample Masisardi ct al. (2011) found, at 20 GHz, a variability 
of ~ 9 per cent at 20 GHz on a timescale of abo ut 9 months, 
slightly larger than what has been found by ISadler et al.l 
(200fi) (6.9 per cent) over a one year timescale for a 100 
mjy fiux density limit sample. Our result is consistent with 
earlier indications (i.e., Giaramella et al. 2004) that the vari- 
ability increases with the time lag, for lags of up to several 
years. 

An inspection of Fig. [7] suggests that the disper- 
sion a„ may be somewhat lower at the fainter fiux den- 
sities. For the fiux density interval 0.2 ^ 5at20G < 
0.4 Jy we find (logio(5'pACo/5'AT20G)) = -0.03 and 
a[logio(5'pACo/5'AT20G)] = 0.17l°;°2. while for 5'aT20G ^ 
0.4 Jy we have (logiQ(S'pAco/5'AT20G)) = -0.05 and 
o-[logio(5'pAco/5'AT20G)] = O.^l^l- The 68 per cent 
confidence errors on a have bee n computed following 
iDanese. de Zotti. fc di TuUiol (|198G| '). The difference is not 
statistically significant, as confirmed by a KS test that yields 
an 3.8 per cent probability that the two distributions come 
from the same parent population. 

We have made the same analysis at lower frequency by 
comparing Spaco^.tghz with the fiux densities at 4.8 GHz 
measured in the AT20G follow-up observations. The least 
square linear relation between the two sets of observations 
(see Fig. El) gives 5pAco[Jy] = (0.92 ± 0.03) • 5'at2og[J2/] + 
(0.05 ± 0.04), confirming the pre- and post-GABB com- 
patibility also at these frequencies. The variability (a^,) is 
38 per cent and the mean and standard deviation of the 



The PACO faint sample 7 



Table 3. Fractions of PACO sources with ~ 19 GHz (more precisely, at 18.768 GHz) flux densities (mjy) within the bins indicated in 
the first column, found in the 32.768 (upper panel) and 39.768 GHz (lower panel) bins listed in the first row. The values in parenthesis 
are the number of sources in each 19 GHz bin. 



•S'lQGHz 


200-250 


250-400 


400-750 


750-1000 


1000-1500 


1500-2000 


2000-5000 


200-250 (23) 


0.11 


0. 


0. 


0. 


0. 


0. 


0. 


250-400 (38) 


0.45 


0.39 


0. 


0. 


0. 


0. 


0. 


400-750 (68) 


0.03 
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Table 4. Source counts at 32.744 (second column) and 39.768 
GHz(third column), with their errors, obtained by extrapolating 
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source counts. 
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Figure 9. Distribution of logjQ(S'pACo/'S'AT20G) at ~ 20GHz. 
The dotted line is the Gaussian best-fit. The dashed lines are the 
mean and the mean ±(t of the distribution. 



logig(S'pACo/SAT20G) distribution are -0.006 and 0.15 re- 
spectively. 



4.3 Source Counts 

We use the simultaneous PACO measurements to estimate 
the distribution of flux densities at 33 or 40 GHz of sources 
with 19 GHz flux density within a given interval. To increase 
the statistics at Sat20G ^ 500 mJy, we have taken into ac- 
count also the PACO bright sample that covers an area a 
factor of 4 larger than the PACO faint sample. Using the re- 
sults in Table [3} we d erive from the 20 GHz AT20G counts 
l|Massardi et alll2010l ). the counts at 33 and 40 GH z. The 
errors have been computed following iGehreld l|l986h . 

Table [3] suggests that additional, albeit small, contribu- 



tions to the counts in the lowest and in the highest 33 or 
40 GHz flux density bins could come from sources outside 
the flux density range covered by the PACO samples. Thus, 
strictly speaking, the derived counts in these bins are lower 
limits. While at low flux densities the incompleteness arises 
from the flux density limit of the PACO faint sample, at high 
flux densities it may be due to the limited sky coverage. 

Figure [TD] compares the present estimates of source 
counts (listed i n Table [H) with those obta ined from the 
Planck ERCSC jPllanck Collaboration|[2011il 'l and from the 
DASI experiment ( Kovac et al.ll2002l ). The agreement is very 
good in the region of overlap, and our data extend the counts 
downwards in flux d ensity by a substantial factor compared 
to those obtain ed bvlPlanck Collaboration! (|2011il ). The pre- 
dictions of the Ide Zotti et al.l 1 20051 ) model are nicely con- 
sistent with these data. 

The extension of the counts downwards in flux density 
is important to quantify the fluctuations in CMB maps due 
to sources below the detection limit. Si. Since the differential 
counts below the completeness limits of the ERCSC (~ 1 Jy 
at 30 GHz and ~ 1.5 Jy at 44 GHz, iPlanck Collaboration! 
(j2011il )) go roughly as dN/dS oc S~^'^, the power spectrum 
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Figure 10. Differential source coun ts at 33 and 40 GHz (squares) 
compared with tlie Planck counts l lPIanck Collaboration! boilili 
at 30 and 44 GH z (asterisks), and with the DASI estimate 
llKovac et al.|[20o3) at 31 GHz (region bou nded by dotted lines). 
The solid lines show the predictions of the lde Zotti et al 



power spectrum at 33 GHz 
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Figure 1 1 . Contribution of the extragalactic point sources to the 
CMB power spectrum (black solid line) as a function of multipole 
(1) for different flux density limits. In each panel the upper (dot- 
dashed) line refer to the ERCSC completeness limits (1 Jy at 33 
GHz and 1.5 Jy 40 GHz), while the lower (dashed line) refers to 
the flux density limit of the PACO faint sample (200 mjy). 



of Poisson flu ctuations due to sources fainter than 5"; goes 
as (|de Zotti e t al. 1996): 



C = 0. 



— — b db oc bi , 
db 



(3) 



where Q. is the beam solid angle. This means that ^ 70 
per cent of the residual point source fluctuations in Planck 
30 and 44 GHz maps, after subtraction of directly detected 
sources, are contributed by sources in the flux density range 
covered by the PACO survey. Note that fluctuations due to 
clustering, not taken into account in the above equation, are 
negligible in the case of radio sources because of the strong 
dilution due to the broadness of the luminosity function. 
The effect of an assessment of the counts down to fainter 



flux limits on the contribution to the power spectrum of ra- 
dio sources below the detection limit is shown in Fig. 1111 
Here the dotted line refe rs to the com pleteness limits of the 
ERCSC (Planck Collaboration! bOllcl). 1 Jy a t 33 GHz and 
1.5 Jy at 40 GHz (|Planck Collaboration!!2011i! ). The present 
assessment of the counts down to 200 mJy decreases the am- 
plitude of the extragalactic source power spectrum as shown 
by the dashed lines. 



5 CONCLUSIONS 

The PACO faint sample is a complete sample comprising 
159 AT20G sources with 5'2oghz > 200 mJy. The sources 
have been observed between July 2009 and August 2010 
in 3 frequency bands of 2 x 2 GHz each, centered at 5.5, 
9, 18, 24, 33, 39 GHz. Most observations were carried 
out with in 10 days from Planck obs ervations. The Planck 
ERCSC jPlanck Collaboratioril!201ld ) reports the detection 
at least at one frequency of 57 PACO faint sources. Of these, 
25 have measurements at more than 5 Planck frequencies. 
Some examples of spectra obtained combining PACO and 
PZancfc observations are shown in Fig.|3l There are 4 ERCSC 
sources within the PACO faint area not included in our sam- 
ple. One is Fornax A, missed by the AT20G because its 
emission extends over an area larger than the ATCA field-of- 
view. The other 3 are at low Galactic latitude and were not 
detected by the AT20G survey. They may be either extended 
Galactic sources or spurious (peaks of the Galactic diffuse 
emission). Thus the fraction of potentially spurious ERCSC 
sources in the PACO faint area is 5%. Many more sources 
should be detected by Planck as new surveys are completed 
and the data analysis improves. Additional constraints on 
the source spectra at Planck frequencies will be extracted 
from Planck maps, when they will become available, by a 
stacking analysis. 

A comparison with AT20G measurements, carried out, 
on average, a few years earlier, has demonstrated that, on 
these timescales, our sources show a rather high variability 
with an rms amplitude of ~ 40 per cent at 20 GHz. 

The analysis of the spectra, within the PACO frequency 
range, of the 98 sources in the PACO faint sample with 
^oGHz < 500 mJy (i.e. below the flux density limit of the 
bright PACO samp le) shows no significa nt differences with 
the bright sample (iMassardi et al.l!201l! ). The main differ- 
ence is an increase of the percentage of steep sources from 
3.6 per cent in the bright sample to 13.3 per cent in the faint 
one. The high-frequency steepening found bv !Massardi et al] 
(2008) and by Massardi ct al. (2011) is confirmed. 

Our data have allowed us to extrapolate the 20 GHz 
source counts to 33 and 40 GHz, down to ~ 200 mJy, i.e. 
a factor > 5 below current estimates from Planck data at 
nearby frequencies (30 and 44 GHz) . Our counts are in very 
good agreement with the Planck ones in the region of over - 
lap and are well accounted for by the !de Zotti et al] (!2005! ) 
model. The assessment of the counts down to faint flux den- 
sity limits substantially improves the control of fluctuations 
due to unresolved sources in Planck maps. 
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